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Abstract 

We determine the transfer functions of two kinds of filters that can be used 
in the detection of continuous gravitational radiation. The first one optimizes 
the signal-to-noise ratio, and the second reproduces the wave with minimum 
error. We analyse the behaviour of these filters in connection with actual 
detection schemes. 
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I. INTRODUCTION 



The detection of gravitational waves (g.w.) is one of the most fascinating and challenging 
subjects in Physics research nowadays. Besides checking the General Relativity theory, the 
detection of this phenomenon will mark the beginning of a new phase in the comprehension 
of astrophysical phenomena by the use of gravitational wave astronomy. 

Although these waves were predicted at the beginning of the century |T|], the research on 
their detection only started around 1960, with the studies of Joseph Weber [Q. The major 
obstacle to this detection is the tiny amplitude the g.w. have [Q. Even though the more 
sensitive detector now operating |^ is capable to detect amplitudes near /i ~ 5 x 10^^^, 
this value must be decreased by several orders of magnitude so that impulsive waves can be 
detected regularly. 

On the other hand, the discovery of pulsars with periods lying in the milliseconds range 
stimulated the investigations on the detection of gravitational waves of periodic origin. Al- 
though these waves generally have amplitudes even smaller than those emitted by impulsive 
sources, periodic sources are continuously emitting gravitational waves in space and they 
can be detected as soon as the correct sensitivity is reached. Since many of the resonant 
mass antennae now operating are designed to detect frequencies near 1000 Hz, the millisec- 
ond pulsars will probably be detected if these antennae ever become sensitive to amplitudes 
h < 10^^''. This value is bigger if we consider the Crab pulsar (/ ~ 60Hz): h ~ 10~^^. 
There is a resonant mass detector with a torsional type antenna (CRAB IV) being developed 
by the Tokyo group P to detect gravitational waves emitted by the Crab pulsar. This group 
expects to reach h ~ 10~^^ soon. 

The main purpose of this paper is a contribution towards the increase in sensitivity of 
resonant mass continuous gravitational wave detectors looking at the use of adequate filters. 
We study two kinds of filters, the first optimizes the signal-to-noise ratio (SNR), and is 
normally used in the detection of impulsive waves 0. The second filter reproduces the 
wave with minimum error. Both filters apparently were not investigated in the continuous 
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gravitational wave context yet. 



II. THE FILTER THAT OPTIMIZES SNR 

Linear, stationary filters obey the relation 

/oo 
k{t')I{t - t')dt'. 
'OO 

k{t) is the impulse response function that characterizes the filter /C , I{t) is the input at the 
filter and 0{t) is the filter output. 

Generally^ X(t) has a useful part, f/(t), and an unwanted part, N{t): I{t) = U{t) + N{t). 
We have a similar relation for the filter output, given by 0{t) = U'{t) + N'{t). 

It is well known from noise theory that the filter /Co that optimizes SNR at its output^, 

must have the following transfer function: 

i^oH = e--"-^, (2) 

with 

/oo 
e-"^'k{t)dt. 
-oo 

to is the instant in which the observation takes place, U{uj) is the Fourier transform of U{t) 
(* denotes complex conjugation) and Sn{u) is the noise power spectrum density: 

/oo 
e"'*"* < N{t)N{t - r) > dr. 
-oo 

The maximum SNR at the optimal filter output is given by the expression 



^We are only considering random, stationary processes in this work. 
^< f{t) > represents the average value of f{t). 
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/Iff-' 

From (0) and (|^) we conclude that a very weak signal will leave the filter when the noise 
is much stronger than the useful signal at the relevant frequency range. 



III. QUASI-MONOCHROMATIC SIGNALS 

Equation (H) is valid as long as U{uj) is well behaved. For example, if U{t) were a strictly 
monochromatic wave like 

U(t) = /lo cos Wot, (4) 

it would be difficult to build this filter since U{uj) = h^diuj — ujq). 

In order to use the optimal filter in continuous gravitational wave detectors we will 
describe these waves as quasi-monochromatic useful signals. It means that the waves that 
reach the antenna will be of the formQ 

2/ioe-"*cos(cJot + t) ,t > 
I cos(cjot + f ) , t < 

The constant a is related to the signal spectral density bandwidth, Au;^ = and the 
corresponding spectral density is of the formQ 



/?2 



2a 2a 
+ 



+ (cj - tuo)^ a"^ + {oo + 000^ _ 
The signal (^) is quasi- monochromatic whenever a << ujq, ujq being its central frequency. 
Note that when a ^ we recover (|^), the monochromatic case. 

The continuous gravitational waves emitted by periodic sources can be regarded as quasi- 
monochromatic waves. The frequency of the Crab pulsar, for example, which is centered 



^We suppose, for simplicity, that the g.w. has only one polarization, namely "+". 
^At this time we admit uj real and luq > 0. 
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near QOHz, has a slow down rate of ~ O.OlHz/year. Besides, the orbital motion of the 
Earth causes a maximum variation of ±0.03Hz/year, and the spinning motion of Earth 
implies a maximum variation of ±2 x 10~^ Hz/ day P). 

For future use in the optimal filter expression, (^, we write the Fourier transform of the 
quasi- monochromatic signal (^: 

l+i 1 



h{uj) = \/2ahQ 



+ 



a'^ + {uj - ujoy + (cu + cuq) 



(6) 



IV. THE MATHEMATICAL MODEL OF THE DETECTOR 

A resonant mass detector can be represented by the scheme of figure |l[ In this model 
F{t) represents the gravitational interaction force between the g.w. and the antenna. The 
two-port circuit is related to the massive antenna and the transducer^, and it is described 
by its admittance matrix yij{uj), which relates the force /i and the velocity vi at the input 
port to the current / and the voltage V at the output port: 

i{iu) = yiifi{uj) + yi2Vi{uj) (7) 

V{uj) = y2ifiiuj) + y22Vl{uj) 

The transducer and the amplifier have force and velocity noise generators represented 
by the stochastic, stationary functions f{t) and v{t), respectively. Sfiuj) [S^ioj)] is the 
spectral density of /(t) [v{t)\. We will assume that these functions are not correlated, so 
that Sf^ioj) = 0. 



^We will adopt a non-resonant transducer Q. 
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FIGURES 



FIG. 1. Resonant mass detector scheme. 

In this model the optimal filter follows the lock-in amplifier. In figure the elements 
that precede the optimal filter in the detector are redrawn ||10|| . 



FIG. 2. Model to the antenna, the transducer and the amplifier. 

In this figure m is the antenna effective mass, k its elastic constant and ^ the damping 
constant. represents the mechanical dissipation at the antenna, f{t) is associated to 
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the transducer back- reaction on the antenna and v{t) represents the wideband, serial noise 
introduced by the amphfier. 

The equation of motion of the system in given by 

mz{t) + fiz{t) + kz{t) = + f{t)] - F{t), 

where z{t) represents the displacement. However, in the calculations we will deal with the 
velocity z(t); at the amplifier output we have z'{t) = z(t) + v{t). 

In the absence of F{t) we can obtain the velocity noise spectral density: 



At the denominator of this expression, the minus signal is used when u; > 0, and the plus 
signal is used when < 0. 

Since the thermal and the back-reaction noises are white noises, the force spectral den- 



sities generated by them obey the following Nyquist relations p!T[: SAuj) = 4 '"^'^'^'^ and 
Sf{uj) = 4 "^^^g^/ _ these expressions, and Tf are the time constants of mechanical loss 
at the antenna and of electrical loss at the transducer and amplifier, respectively. They 
are related to the energy decay time of the antenna, = according to the expression 
= ^ + ^1 where Qq is the antenna quality factor. is the antenna temperature and Tf 
is the back-reaction noise temperature, ks is the Boltzmann constant. 

The function Ss{uj) that appears in (P) represents a serial white noise introduced in 
the useful signal by the electrical network, and it has the following expression: Ss{i^) = 
j^j^ARtkBTr, where ?/22 comes from (J^). Rt is the real part of the impedance at the trans- 
ducer output and T,. is the circuit noise temperature. We assume that the amplifier has a 
large but limited bandwidth, as it occurs in practice , given by r^^. Generally 

T-o » Ta. (9) 

Introducing the antenna's equivalent temperature, Tg, defined by the relation ^ = ^ + ^, 
equation (§) becomes 



S'Ar(t^) = — — y + , (10) 

This is the complete expression for the total noise spectral density at the filter input. 

The velocity that the g.w. generates at the antenna in the absence of the noises 
and f{t) has the following Fourier transform: 

2m 1 + tTo[uj — ooq) 

F(uj) is the Fourier transform of the g.w. force on the antenna, which is given by the relation 

m 



nt) = E qap^.h^p. (12) 

a,/3=l 

Qais is the dynamic part of the mass quadrupole tensor of the antenna. It is a matrix of 
constant elements which depend on the antenna's geometry and mass distribution. For our 
calculations we use Qxx ~ PiiY: where / is the characteristic length and p is the density of 
the antenna. 

Equation (pA]) is the useful signal contribution at the filter input, which corresponds to 
the spectral density]^ 

\2mJ l + ro(a;-u;o) 
Sf{uj) is the spectral density of the force (0) and it has the form 

Sf{u;) = \qlu;',S,{u;). (14) 



V. THE FILTER THAT OPTIMIZES SNR FOR THE CONTINUOUS G.W. 

DETECTOR 

Using (^, (^) and ( [TT| ) and adopting to = we obtain the transfer function of the filter 
that optimizes SNR for the model considered in the preceding section: 



^We will simplify the expressions adopting to > henceforth. 
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KJio) 



+ 



(15) 



To simplify this expression we introduced the following definitions: Q = ^^^^f^, Np = 
2TokBTjm and A^^ = 4:RtkBTr/\y22\'^- 

The maximum SNR related to this filter is obtainable from (^, (ID) and ([TB|), and it 
corresponds to 



1 poo 

SNR„ = — 



(t 



Note that when a « Tq-^ << r^^, this expression becomes 



■duo. 



(16) 



Assuming, for simplicity, that Ng = Nj^, we find 



(17) 



6AkB a niTe 

By imposing SNRo > 1 we obtain the following condition on the parameters of the detector: 



> 8.83 X 10"^^ 



(18) 



To illustrate the use of the filter ( [T^ ) we will adopt two different realistic detectors. 
In both cases we will assume that A^^ = Np and tq = IO^Tq (see equation (^. The first 
of them, detector A, is designed to detect Crab pulsar {ujq = 376.99Hz, Hq ~ 10"^'' and 
a = 3.17 X lO^^^Hz. This bandwidth arises from the slow down of the pulsar after 100 
milliseconds of observation); this detector has the following characteristics: Qq = 5 x 10'', 
m = 1200% §, Te = 5K, p = 2.7Ag.cm-^, I = 1.1m. 

Under these conditions, (|1^) shows a very narrow peak centered in uq. It implies a 
maximum SNR at its output given by SNRo = 1. If the signal bandwidth is smaller 
(implying a smaller observation time), we have SNRo > 1. 



^We are supposing that the back-reaction noise and the serial noise give the same contribution to 
the total signal at the filter input. 
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On the other hand, if detector A has a lower equivalent temperature we attain SNRo = 1 
with a longer observation time. For instance, we obtain this result if Te ~ O.OSi^', and 
a = 3.17 X 10^^ Hz. Such bandwidth corresponds to a 10 seconds observation time of the 
Crab pulsar's slow down . 

The other detector considered, detector B, is designed to detect the millisecond pulsar 
PSR 1937+214 H (too = 8,066A7Hz, ho ~ lO^^^). Since we did not find any information 
about the bandwidth of the g.w. emitted by this pulsar, we will assume that it has the value 
a = 1.2 X lO^^Hz. This detector is characterized by Qo = 5 x 10^, p = 2.74g.cm~^, I = 2m, 
m = 2300kg and Te ~ O.li^'; these are typical values of several ultracryogenic cylindrical 
antennae. 

For detector B, (|l^) also shows a very narrow peak centered in uq, implying SNRo = 1. 
Like detector A, SNRo is greater if the signal bandwidth is smaller. 

VI. THE FILTER THAT REPRODUCES THE USEFUL SIGNAL WITH 

MINIMUM ERROR 

After the g.w. is detected we have to determine its shape with minimum error. This can 
be accomplished with the help of an adequate linear filter, /C^, designed to reproduce the 
useful signal with the greatest possible accuracy (depending on the noise and the useful signal 
present at its input). This accuracy is characterized by the mean square error, < e^(t) >, 
which is obtained from the instantaneous reproduction error, e{t), defined by 

e{t) = Oit)-v{t). (19) 

ri{t) is the desired signal at the filter output. In a simple filtering process, as the one we 
are considering in this work, rjit) must be equal to the useful signal at the filter input, U{t). 

We obtain the transfer function Kr{uj) of the filter /C^ by imposing < e^(t) >=< 
e^(^) >min- This condition implies [0 
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supposing there is no cross-correlation between U{t) and N{t). The corresponding mean 
square error is 

9/ N 1 Su(uj)SNiuj) , 

..^<.'W ^X. swmUm ''" '''' 

From this equation it is evident that the error becomes smaller if so becomes the noise. 
On the other hand, if the noise is too strong {Sn{^^) oo) it results Kr{uj) — >• and no 
signal leaves the filter. 

If we define the total power of a signal Z{t) by < Z'^{t) >= ^ Sz{uj)duj, the SNR at 
the filter input will be0 

SNR,nput = T^^y^- (22) 

At the filter output the signal Z{t) will have the following total power: 

1 1"°° 

<Z'^{t)>=— \Kr{uj)\^Sz{uj)duj. (23) 

ZTT J -oo 

Using this relation we can find the SNR at the filter output, 

SNRoutput = < ^/2jij > • (24) 

VII. SIMPLE FILTERING OF THE QUASI-MONOCHROMATIC SIGNAL 



We now consider the particular case of noise spectral density given by ([10| ) and useful 
signal spectral density given by (|lT]). In this case the filter that reproduces U{t) with 
minimum error has the following transfer function: 

K,.{u) = I 1 + }±^E^^K_}±tz^f^ I . (25) 

[l+a-2(c^-a;Q)2][l+r^(<^-a;o)2] 



*Note that ( p2[ ) is different from (|T]). This happens because we are now interested on the total 
spectrum of the useful signal, while in the analysis of the first kind of filter we were interested only 
on the maximum amplitude of this signal. 

11 



The minimum error introduced by this filter is obtained from using (0) and ([131): 

\ l + (w-c^o)V|.f ^ l+{(^-a;o)V2 / 

In detector A the total noise power at the filter input is < N'^{t) >= 5.8 x lO^^^m^/s^. 
If we use filter ( pS] ) in this detector we obtain = 3.42 x which is ~ 6 x 10~^ 

times smaller than < N'^{t) >. Comparing (|2^ ) and for this case, we find SNRmput ~ 
8.34 X 10'« and SNRautput ~ 0.68, so that SNRoutput ~ 8.11 x lO^SA^i^i^p^f. 



On the other hand, using filter (^Sf) in detector B we obtain an output error of er = 
6.24 X 10~^^m^/s^, which is almost 10~^ times smaller then the input noise, < N'^{t) >= 
6.06 X 10-26mVs2. In this case, SNRinput ~ 1-46 x 10"^ and SNRoutput ~ 0.67, which imply 
SNRoutput ~ 4.67 X 10^ SNR^nput. 

If the Crab bandwidth were a ~ 1.85 x lO^^^Hz we would obtain SNRoutput ~ 1.05; 
this bandwidth also implies er ~ 2.94 x 10~'^^m^/s^. We would find the same SNRoutput 
if the bandwidth of PSR1937+214 were a ~ 6.9 x 10~^Hz, corresponding to ~ 5.36 x 

10~3W/s2. 



VIII. CONCLUSIONS 

We have derived expressions for the transfer functions of two kinds of filters, both de- 
signed to detect continuous monochromatic waves. The first filter, JCo, optimizes SNR at its 
output (equation (|l])) and is important for a first detection of the wave. The second filter, 
JCr, reproduces the wave with minimum error and should be used when we intend to know 
the complete shape of the wave. 

In the study of JCo we have first analysed the detection of Crab pulsar. Supposing 
Np = Ng (see section V) we concluded that this pulsar could be detected if its signal were as 
monochromatic as ~ 3 x lO^^^Hz; it means a sampling time of the order of 100 msec. We 
have also analysed a possible detection of PSR 1937+214 and suggested a maximum limit 
for its bandwidth, allowing its detection by third generation resonant mass detectors. 
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The use of /Cj. is effective in the detection of both pulsars if their signals are still more 
monochromatic. We then obtain SNRoutput > 1 with the form of the signal as preserved as 
possible. 

With the continuous optimization of present detectors the condition of monochromaticity 
of the signal becomes weaker. For example, if Np > Ng, the bandwidths of the signal can 
be larger than those we obtained in this paper. Besides, the inequality (|18|) can be used 
as a reference to optimize resonant mass continuous gravitational detectors. Note that 
continuous sources with high frequency, small bandwidth and high amplitude are the most 
favourable for detection with less improvement of the detector. On the other hand, the 
detector should have a small equivalent temperature and materials with high quality factor 
and density should be preferred for the antenna body, which should be as large as possible. 
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